Elementary Charge Transfer Processes in a Superconductor-Ferromagnet Entangler 
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We study the production of spatially separated entangled electrons in ferromagnetic leads from 
Cooper pairs in a superconducting lead. We give a complete description of the elementary charge 
transfer processes, i) transfer of Cooper pairs out of the superconductor by Andreev reflection 
and ii) distribution of the entangled quasiparticles among the ferromagnetic leads, in terms of 
their statistics. The probabilities that entangled electrons flow into spatially separated leads are 
completely determined by experimentally measurable conductances and polarizations. Finally, we 
investigate how currents, noise and cross correlations are affected by transport of entangled electrons. 
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A solid state entangler is an eleetronie analog of the op- 
tical setups used for experimental Bell inequality tests, 
quantum cryptography and quantum teleportation [ij. 
Ideally, such a device should produce separated currents 
of entangled electrons. Superconductors are suitable can- 
didates as sources in solid state entanglers since Cooper 
pairs constitute entangled states. This prospect has mo- 
tivated several papers addressing the properties of hybrid 
superconductor and normal metal entanglers 0, 0, |3j 01 ■ 
One of the challenges is to prevent processes where pairs 
of entangled particles reach the same lead, i.e. are not 
spatially separated. Electrons from Cooper pairs are en- 
tangled in spin and energy space, and separation of pairs 
into different leads using ferromagnets or quantum dots 
has been suggested Upon filtering, only the spin or 
energy part of the two-particle wave function collapses, 
depending on whether ferromagnets or quantum dots are 
used. Respectively, energy or spin entanglement remains 
Here we consider separation by ferromagnets. 

Solid state entanglers have been analyzed in Refs. 
[3, H, 0, [l] in terms of currents, noise and cross correla- 
tions. A more direct approach, describing the elementary 
charge transfer processes in terms of experimentally con- 
trollable parameters is certainly desirable. We demon- 
strate how this is possible through the full distribution 
of current fluctuations, the full counting statistics (FCS), 
of the solid state entangler 0, 0, Q ■ The FCS provides 
complete information about currents, noise, cross cor- 
relations and higher cumulants, and even more impor- 
tantly, allows direct access to the probability for transfer 
of charge between different parts of the device. 

We consider the singlet superconductor-ferromagnet 
(S-F) device shown in Fig. [T] A normal metal cavity (c) 
is connected to one superconducting terminal and several 
ferromagnetic terminals via tunnel junctions. The cavity 
is under the influence of proximity effect. In this device, 
charge transport occurs via two processes: i) Transfer 
of Cooper pairs out of the superconductor by Andreev 
reflection and ii) distribution of the entangled quasi- 




FIG. 1: Circuit theory representation of generic F-S entangler. 
Entangled electrons from the singlet superconductor S enter 
the cavity c through a tunnel barrier with conductance gs and 
escape through ferromagnetic interfaces with conductance <?„ 
and spin polarization \gn\/gn into drains F„. Arrows indicate 
magnetization directions gji/lgnl- 



particles among the ferromagnetic leads. The distribu- 
tion can occur via Direct Andreev (DA) reflection, where 
a entangled pair is transferred into lead F„ or crossed 
Andreev reflection (CA), where each particle of the en- 
tangled pair is transferred into spatially separated leads 
Fjn and F„ (n ^ m). CA produces spatially separated 
entangled electrons. Since the ferromagnetic terminals 
are at the same voltage and we consider zero tempera- 
ture, there is no direct electron transport between the 
ferromagnetic terminals Q . 

Our general results for the counting statistics show 
that the processes i) and ii) are independent and therefore 
the statistics can be factorized. This novel factorization 
and the probability distribution for process ii) reveals 
the precise dependence of the probabilities for CA and 
DA processes on the experimentally measurable conduc- 
tances and polarizations of the ferromagnetic leads. We 
flnd that the probability that the electrons of a Cooper 
pair transferred into the cavity are detected in terminals 
m and n is 



(1) 



where g = X]„ 5n and g = J^n The probabilities Pmn 
depend solely on the conductances g„ and spin polariza- 
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tion conductances g„ of the ferromagnetic leads which 
can be determined by magnetoresistance measurement 
in the normal state. Eq. ^ shows that the detected 
two-particle processes originate from a pure spin singlet 
density matrix subspace [10|]. We emphasize that ((T|) in 
combination with the Cooper pair transfer probability, 
to be discussed below, allows for an unambiguous identi- 
fication of all statistical properties of the charge transfer. 
Using the magnetization dependence of the probabilities 
one can violate Bell's inequality straightforwardly and 
demonstrate entanglement. In a device with e.g. two 
ferromagnetic terminals, the probability to separate the 
entangled quasi-particle pair into different leads is en- 
hanced in an antiparallel magnetization configuration. 

It is quite remarkable that all statistics (noise, cross 
correlations, and higher order cumulants) in this device 
are completely determined by the normal state transport 
conductances between the cavity and the terminals and 
no additional parameters need to be introduced. The 
relative orientations of the magnetizations are the control 
parameters in an experimental situation. The fraction of 
the CA current and, therefore, the spatially separated 
entangled pair currents follow from these. 

Ferromagnet detection of entangled spin singlets from 
a ballistic normal conductor was considered in Rcf. [lo| . 
In that device, there are also one-particle transfers, which 
can contribute substantially to the current and the noise. 
Another important qualitative difference between our de- 
vice and the system in Ref. is the latter's strong cou- 
pling to the detectors which can distort the spin singlets 
emitted from the source and induce triplet correlations 
upon detection. Also, in the limit of a weak coupling to 
the detectors, there are no two-particle processes in the 
system of Rcf. [l^l- In contrast, the electrons in our S-F 
device are always detected from a pure spin singlet state. 

The charge transfer probabilities arc obtained by iden- 
tifying the elementary processes in the many-body charge 
counting statistics. The statistics is determined by the 
cumulant generating function (CGF) S{xi,...,Xn) = 
S{{xn}) of the probabihty P {{Nn}) to transfer in a time 
interval to, Ni electrons to Fi, N2 electrons to F2, and 
so on. Our main finding is the statistics 
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for Nn even and positive. The interpretation of this 
result is that the charge transfer is given by two indepen- 
dent processes. The first factor Ps{2N) is the probability 
that N = X]n^n/2 Cooper pairs are emitted from the 
superconducting source terminal into any of the detec- 
tors. The second factor P ({iV„}|2A^) in (l2b| is the con- 
ditional probability that Nn out of the 2A^ electrons have 
been transferred into ferromagnetic terminal F„. Below 
we will explain in detail how our calculation yields con- 



crete expressions for the elementary processes described 
by P {{Nn}\2N). These results facilitates a unique inter- 
pretation of the transfer of spin singlet electron pairs. 

To complete the full statistical description, we now 
supply the microscopic expressions for the two probabil- 
ities in (|2bp . The Cooper pair transfer probability is ob- 
tained from Ps{2N) = / dx/(27r)exp(S's(x) - iNx) with 
a CGF Ssix) given by 
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g^ + g^. The contact to the supercon- 



ducting terminal is characterized by a spin-independent 
conductance gg. The 7r-periodicity of >S'g(x) on x en- 
sures that an even number of charges is transferred. The 
2A^ electrons are distributed among the F„ terminals ac- 
cording to the multinomial distribution P {{Nn}\2N) = 
J d"x/(27r)"exp(5N({Xn}) - X«iV«) with a CGF 
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The concrete form of the two-particle probabilities pmn 
to detect one charge in terminal m and one in terminal 
n is given in ([1]). We will explain below how our calcula- 
tion based on the circuit theory of full counting statistics 
determines Eqs. ([!]), ([31), and g]). 

We emphasize that our interpretation in terms of two 
independent processes of charge transfer is based on a de- 
tailed calculation with a general result for the counting 
statistics. We have not made any a priori assumptions 
on the initial state of the superconducting source termi- 
nal or the ferromagnetic terminals, except that they are 
reservoirs at zero temperature with a voltage bias eV 
applied. The direct result of our calculation is the CGF 
S{{xn}) of the S-F entangler of Fig. [1] Its explicit ex- 
pression is found in ([3]) by replacing the factor e^*-^ with 
exp{S'N({xn})/A^} given in (g]). The factorization in ^ 
can be proven straightforwardly from S{{xn})- Actually, 
such a factorization is valid for any CGF where the x de- 
pendence is exp{S^{{xn}) / N} , irrespectively of its form 
or the probabilities Pmn- 

We now discuss some consequences of the charge count- 
ing statistics. FCS enables us to express the current 
and noise correlations in a compact and meaningful form. 
The currents /„ = iie/to)dS({xn})/dxn\x„.=o are 



-, In ^ Ipn, (5) 



where p„ = l^mPmn is the probability to detect one 
of the electrons in terminal n, irrespective where the 
second electron goes. The combined probabilities can 
be directly accessed in the noise correlators between 
current fiuctuations in terminals m and n, Cmn = 
(-2eVto)a25({x4)/5Xm%n|x„„„=o: 
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2(1 - F2)pmPn] . (6) 
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The Fano factor for Cooper pair transport is defined as 
the ratio of the full current noise C = X^mn ^mn to the 
Poissonian noise of doubled charges, = C/4eI, and is 
explicitly found 2(1 - F2) = [5 - g^/ig^ + g^)]x/il + 
where x = <?g/(5^ — g^)- These expressions for the current 
and the noise provide a transparent interpretation of the 
transport processes. The current in ([5]) is proportional to 
gsig^—g^), since two particles have to tunnel through the 
double junction to transfer a Cooper pair from S. The de- 
nominator is due to the proximity effect 11 , l3, l3| and 
enhances the current drastically in comparison to cal- 
culations based on the tunneling Hamiltonian The 
current into each terminal is then weighted accord- 
ing to the probability Pn- We might also distinguish the 
contributions to the current originating from crossed and 
direct Andreev reflection. The probability to detect a DA 
reflection in terminal n is given by Pnm and the proba- 
bility for CA detection in different terminals to 7^ 71 is 
given by Pmn ■ We find the ratio of the crossed current to 
the total current as 



-^n Pn Pnn 
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This ratio is independent of the coupling to the super- 
conducting terminal. We further observe that the crossed 
current is enhanced by increasing the polarization of the 
contact n and is additionally favored by aligning the 
magnetization g„ opposite to the average magnetization 
g. These results are a direct consequence of the spin- 
singlet nature of the Cooper pairs. Enhancing the mag- 
nitude of the polarization Ignl/gn of one terminal reduces 
the total current, but enhances the crossed part of the 
Andreev current, since the tunneling of one spin-singlet 
electron-hole pair through the same contact is strongly 
suppressed. 

The sign of cross correlations in three-terminal beam 
splitters has been considered for various devices both ex- 
perimentally [3] and theoretically 0, [l3|. Studies of 
noise [lB| and FCS [l3| for a beam splitter with entan- 
gled electrons show that entanglement gives qualitatively 
different noise characteristics compared to transport of 
non-entangled electrons. The physical origin of positive 
and negative contributions to the cross correlators can in 
our case be understood from the dependence on the two- 
particle probabilities in ([B]). CA reflection leads to posi- 
tive cross correlations since two particles are transferred 
simultaneously into and F„ (bunching behavior) [l6| . 
A negative contribution (anti bunching) that does not de- 
pend on entanglement, is induced by the fermion exclu- 
sion principle: The transfer of one electron-hole pair into 
F„i by DA reflection, prevents the simultaneous trans- 
fer of another pair into F„. However, if the electron- 
hole pair transfers are not temporally correlated (Pois- 
sonian statistics), the exclusion principle does not affect 
the cross correlations. This is the case when there is 
strong asymmetry in the junction conductances g and gs 



(.9 ^ 5s or 5 ^ gs) so that the Fano factor F2 = 1. In 
this limit the negative contribution —2(1 — F2)pmPn in 
(l6|) vanishes. Scattering matrix calculations give similar 
results for Cmn 

The strongly asymmetric case is particularly interest- 
ing since the cross correlations (to ^ n) Cmn = '^^Ipmn, 
are a direct measure of the probability that electrons from 
a Cooper pair are transferred into different terminals. 




FIG. 2: (color online) Left: Contour plots of F12 for antipar- 
allel and parallel (inset) configuration (|gi| = |g2|). Positive 
regions in green. Right: Plots of F12 in antiparallel configura- 
tion as a function of the conductance asymmetry gs/g, legend 
denotes the value of the polarization \gn\/gn- Red and blue 
horizontal lines in left panel correspond to red and blue curves 
in right panel. 

To illustrate our theory, let us now consider the three- 
terminal version of Fig. [T] with the superconducting 
source terminal S and two ferromagnetic drains Fi and 
F2 . The ferromagnetic magnetizations can in this device 
be utilized as filters to produce currents of entangled elec- 
trons in separated leads. Let us consider |gi | = |g2 1 in the 
following and define Fano factors F^n = Cmn/(2e/). The 
autocorrelation noise -^11(22) will be reduced in antipar- 
allel alignment gi = — g2 as compared to a S-N system 
(g„ = 0) due to enhancement of CA. The cross corre- 
lation F12, shown in Fig. [5] can have both positive and 
negative sign depending on the conductance asymmetry 
gs/g and the spin polarization. The positive contribu- 
tion to F12 is proportional to gig2 + (— )|gi| |g2| in the 
parallel (antiparallel) alignment demonstrating how spin 
filtering of entangled pairs enhances (reduces) the corre- 
lation between currents in Fi and F2 with respect to an 
S-N system Q . Note that for sufficiently large spin polar- 
ization, F12 can be positive for the entire range of gs/g 
in the antiparallel alignment (region above blue line in 
left panel of Fig. [J) , whereas it remains always negative 
in the parallel alignment for gs/g — 1 (inset of Fig. 
The change of sign in F12 by switching from antiparallel 
to parallel alignment is due to the enhanced probability 
of CA events, see ([7]). 

We will finally outline the calculation that yields the 
FCS of the considered devices. We utilize the circuit 
theory of mesoscopic superconductivity 0, [3l and rep- 
resent the circuit in terms of terminals, cavities and con- 
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nectors. Terminals are described by equilibrium qua- 
siclassical Green's function matrices Gn determined by 
electrochemical potential and temperature. Our nota- 
tion for matrix subspace is: ~ for spin, ' for Nambu, 
and " for Keldysh. Pauli matrices are denoted Tj. At 
zero temperature we consider < E < eV where the 
Green's functions for all ferromagnetic terminals F„ are 
G = Tsfa + (fi + 1x2) where V is the voltage of the ferro- 
magnetic terminals and E the quasiparticle energy. The 
superconductor S is at zero voltage and has Green's func- 
tion Gs = Ti, where we assume E A, A being the gap 
of S. The terminals are connected to a cavity c which is 
under the influence of proximity effect from S. The cav- 
ity is described by an unknown Green's function Gc, as- 
sumed isotropic due to chaotic or diffusive scattering. We 
assume that c is large enough so that charging effects can 
be neglected, and small enough so that Gc is spatially ho- 
mogeneous. The circuit theory is formulated in terms of 
generalized matrix currents Ij in spin(KiNambu(XiKeldysh 
matrix space and from the matrix current conserva- 
tion Ij = 0. This determines the Greens function 
on the node together with the normalization condition 
G^ = 1. The matrix currents can have arbitrary struc- 
ture, and allow to derive the FCS by introducing the 
counting fields Xn for each terminal according to Q 
Gn{x) = e^^"^^^^/^Ge~^^"^^'^^^'^ . Spin active connectors 
are taken into account by spin dependent transmission 
and reflection amplitudes tJ? ^ and ^ for particles inci- 
dent on the interface n from the cavity side in channel k 
with spin a. The matrix current through a spin active 
tunnel barrier between c and F„ evaluated at the cav- 
ity side is [13, El /„ = [g„G„/2 + {g„Tf3, G„} /4, Gc] . 
Here, gn = 9Q^k a 1*1- o-P ^^"^ tunnel conductance and 
5q = e'^/h the conductance quantum. The magneti- 
zation direction is encoded in the direction of g„, and 
the conductance polarization in that quantization axis is 

|g«l = ffQl]fc(l*fe,Tp - I^I-.J^)- ^® ^^'^^ neglected here 
an additional term related to spin dependent phase shifts 



upon reflection at the interface [20|, |2l| , as these are small 
in some material combinations or can be suppressed by 
a thin, non-magnetic oxide layer [2^ . The matrix cur- 
rent between c and S is /s = gs[Gs,Gc]/2 [l^. We take 
into account the spin structure of matrix currents /„ and 
Green's functions in S-F-systems, and derive the CGF 
in the linear response regime and for eV <^ A, gen- 
eralizing Ref. [H: S* = io/(4e2) / dEj^p^, where 

{Ap} is the set of eigenvalues of the matrix M deflned 
by writing matrix current conservation in the cavity 
Ij = [M, Gc] = 0. The non-trivial spin matrix struc- 
ture of /„ determines the magnetization dependence of 
transport processes in the system. Carrying out this pro- 
cedure yields the FCS for the setup in Fig. 1. 

In conclusion, we have investigated the elementary 
charge transfer processes of a S-F entangler. Charge 
transfers occurs via two statistically independent pro- 



cesses, i) Cooper pairs are transferred out of the super- 
conductor by Andreev reflection and ii) entangled quasi- 
particles are distributed among the different ferromag- 
netic leads. The probabilities for entangled electrons to 
flow into spatially separated leads are completely deter- 
mined by experimentally measurable conductances and 
polarizations. This allows complete knowledge of the 
statistics of charge transfer in the S-F entangler. 
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